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Abstract 

We deal with the general problem of scattering by open-arcs in two-dimensional space. We 
show that this problem can be solved by means of certain second-kind integral equations of the 
form NS[(p] — f, where N and S are first- kind integral operators whose composition gives rise to 
a generalized Calderon formula of the form NS = Jq +K in a weighted, periodized Sobolev space. 
(Here Jq is a continuous and continuously invertible operator and if is a compact operator.) 
The NS formulation provides, for the first time, a second-kind integral equation for the open-arc 
scattering problem with Neumann boundary conditions. Numerical experiments show that, for 
both the Dirichlet and Neumann boundary conditions, our second-kind integral equations have 
spectra that are bounded away from zero and infinity as k — > oo; to the authors' knowledge 
these are the first integral equations for these problems that possess this desirable property. 
This situation is in stark contrast with that arising from the related classical open-surface 
hypersingular and single-layer operators N and S, whose composition NS maps, for example, 
the function <fi = 1 into a function that is not even square integrable. Our proofs rely on three 
main elements: 1) Algebraic manipulations enabled by the presence of integral weights; 2) Use 
of the classical result of continuity of the Cesaro operator; and 3) Explicit characterization of 
the point spectrum of Jq , which, interestingly, can be decomposed into the union of a countable 
set and an open set, both tightly clustered around —i. As shown in a separate contribution, 
the new approach can be used to construct simple spectrally-accurate numerical solvers and, 
when used in conjunction with Krylov-subspace solvers such as GMRES, gives rise to dramatic 
reductions of Krylov-subspace iteration numbers vs. those required by other approaches. 

1 Introduction 

The held of Partial Differential Equations (PDEs) with boundary values prescribed on open surfaces 
has a long and important history, including significant contributions in the theory of diffraction by 
open screens, elasticity problems in solids containing cracks, and fluid flow past plates; solution to 
such problems impact significantly on present day technologies such as wireless transmission, elec- 
tronics and photonics. From a mathematical point of view, besides techniques applicable to simple 
geometries, existing solution methods include special adaptations of finite-element and boundary- 
integral methods that account in some fashion for the singular character of the PDE solutions 
at edges. With much progress in the area over the last sixty years the field remains challenging: 
typically only low-frequency open-surface problems can be treated with any accuracy by previous 
approaches. 

In this paper we focus on the problem of electromagnetic and acoustic scattering by open 
arcs. In particular, we introduce certain first-kind integral operators and S w whose composi- 
tion gives rise, after appropriate change of periodic variables, to a generalized Calderon formula 
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NS = Jq + K — where is a continuous and continuously invertible operator and where K is a 
compact operator — together with associated second-kind open-surface integral equations of the form 
NS[ip] = f. This approach enables, for the first time, treatment of open-arc scattering problems 
with Neumann boundary conditions by means of second kind equations. Further, a wide range of 
numerical experiments [11] indicate that, for both the Dirichlet and Neumann boundary conditions, 
our second-kind integral equations have spectra that are bounded away from zero and infinity as 
k — > oo, and give rise to high accuracies and dramatic reductions of Krylov-subspace iteration 
numbers vs. those required by other approaches. 

Integral equation methods provide manifold advantages over other methodologies: they do not 
suffer from the well known pollution [2] and dispersion [17] errors characteristic of finite element 
and finite difference methods, they automatically enforce the condition of radiation at infinity 
(without use of absorbing boundary conditions), and they lend themselves to (typically iterative) 
acceleration techniques [4, 9, 29] — which can effectively take advantage of the reduced dimension- 
ality arising from boundary-integral equations, even for problems involving very high-frequencies. 
Special difficulties inherent in open-surface boundary-integral formulations arise from the solu- 
tion's edge singularity [21, 23, 32]. Such difficulties have typically been tackled by incorporating the 
singularity explicitly in both Galerkin [32, 33, 35] and Nystrom [1, 18, 26] integral solvers; with one 
exception (introduced in the contributions [1, 18] and discussed below in some detail), in all of these 
cases integral equations of the first kind were used. While providing adequate discretizations of 
the problem, first-kind integral equation can be poorly conditioned and, for high-frequencies, they 
require large numbers of iterations and long computing times when accelerated iterative solvers as 
mentioned above are used. 

(The literature on the singular behavior of open-arc solutions is quite rich and interesting from a 
historical perspective: it includes the early analysis [31], corrections [5,6] to early contributions [3, 
24] , the well-known finite-energy condition introduced in [25] , the integral equation formulation [23] 
and subsequent treatments for integral approaches for these problems, leading to the first regularity 
proof [32] and the comprehensive treatment [21] which establishes, in particular, that for C°° open 
surfaces with C°° edges, the integral equation solution of the Dirichlet (resp. Neumann) open- 
edge problems equals a C°° function times an unbounded (resp. bounded) canonical edge-singular 
function.) 

As mentioned above, iterative solvers based on first-kind integral equations often require large 
numbers of iterations and long computing times. Attempts have been made over the years to 
obtain second-kind open-surface equations and, indeed, second kind equations for open surfaces 
were developed previously by exploiting the diagonal character of the logarithmic single layer, at 
least for the case of the Dirichlet problem for the Laplace equation [1, 18]. Unfortunately, as shown 
in [11], direct generalization of such approaches to high-frequency problems give rise to numbers 
of iterations that can in fact be much larger than those inherent in first-kind formulations. Efforts 
were also made to obtain second kind equations on the basis of the well known Calderon formula. 
The Calderon identity relates the classical single-layer and hypersingular operators S c and N c that 
are typically associated with the Dirichlet and Neumann problems on a closed surface T c : for such 
closed surfaces the Calderon formula reads N C S C = — 1/4 + K c , where K c is a compact operator in 
a suitable Sobolev space. Attempts to extend this idea to open surfaces were pursued in [12,28]. 
As first shown in [28], there indeed exists a related identity for the corresponding single-layer and 
hypersingular operators N and S on an open surface T. As in the closed-surface case, we have 
NS = —1/4 + K; unfortunately, however, a useful functional setting for the operator NS does 
not appear to exist (K is not compact). As shown in Appendix B, for example, the composition 
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NS maps the constant function ip = 1 on an open surface into a function that tends to infinity at 
the boundary of V like l/d, where d denotes the distance to the curve edge (see Appendix B). In 
particular, the formulation NS cannot be placed in the functional framework put forth in [32, 33, 35] 
and embodied by equations (10), (11) and definition 1 below: a function with l/d edge asymptotics 
is not an element of H~2(F). 

In view of the aforementioned regularity result [21] — which, can be fully exploited numerically 
through use of Chebyshev expansions [1, 26] in two-dimensions and appropriate extensions to the 
high-order integration methods [9] in three-dimensions — , and on account of the results of this 
paper, use of the combination N^S^ enables low-iteration-number, second-kind, super-algebraically 
accurate solution of open surface scattering problems — and thus gives rise to a highly efficient 
numerical solver for open surface scattering problems in two and three dimensions; see [10, 11]. 

This paper is organized as follows. After introduction of necessary notations and preliminaries, 
the main result of this paper, Theorem 1, is stated in Section 2. Section 3 contains the main elements 
of the theory leading to a proof of Theorem 1. Necessary uniqueness and regularity results for the 
single-layer and hyper-singular weighted operators under a cosine change of variables, which have 
mostly been known for a number of years (see [30, Ch. 11] and the extensive literature cited therein) 
are presented in Section 4; the inclusion of these results renders our text essentially self contained, 
and it establishes a direct link between our context and that represented by references [32,33,35]. 
Building up on constructions presented in earlier sections, the proof of Theorem 1 is given in 
Section 5. Two appendices complete our contribution: Appendix A presents a version adequate 
to our context of a known expression linking the hypersingular operator and an integro-differential 
operator containing only tangential derivatives; Appendix B, finally, demonstrates that the image 
of the composition NS of the un- weighted operators N and S is not contained in H~z. 



2 Preliminaries 

Throughout this paper T is assumed to be a smooth open arc in two dimensional space. 



2.1 Background 

As is well known [30, 33, 35], the Dirichlet and Neumann boundary- value problems for the Helmholtz 
equation 

( Au + k 2 u = outside T, u\r = f, feH^(F) (Dirichlet) . . 

[Av + k 2 v = outside T, f^|r = 5, geH~^(F) (Neumann) 

admit unique radiating solutions u,v € H^ oc (M. 2 \F) which can be expressed in terms of single- and 
double-layer potentials, respectively: 



u(r) = y r G fe (r,r>(r')^' (2) 

and 

for r outside F. Here n r / is a unit vector normal to F at the point r' G F (we assume, as we may, 
that n r / is a smooth function of r' E f), and, letting Hq denote the Hankel function, 



GJrr')-i i H o(k\r-r'\), k>0 



(4) 
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and 

dG k (r,r') 

— ^- — = n r / • V r /G fc (r,r). (5) 
Denoting by S and N the single-layer and hypersingular operators 

SM(r)=|G fc (r,r>(ry , r G T, (6) 



and 



def .. d f dG k (r + zn r ,r') , , 

Inn — / ^— -u(r )di , r G T, 

Jr dn r > 



(7) 



z^o+ dz 

the densities fi and v are the unique solutions of the first kind integral equations 

SM = / (8) 

and 

N[u] = g. (9) 
As shown in [32, 33, 35], the operators S and N define bounded and continuously invertible mappings 

S : H~i(r) -»• fl"3(r), and (10) 

N : #3(r) -> fT - 3(r), (11) 
where for seK, the space H S (T) is defined below. 

Definition 1 Let G± be a domain in the plane, with a smooth boundary G\, let seR, and assume 
G\ contains the smooth open curve V . The Sobolev space H S (T) is defined as the set of all elements 
f G H s (Gi) satisfying supp(/) C F. . 

The mapping results (10), (11) provide an extension of classical closed-surfaces results: for a 
closed Lipschitz surface T c and for any s G R, the closed-surface single-layer and hypersingular 
operators S c and N c define bounded mappings 

S c : ^ s (r c ) -»• H S+1 (T C ), (12) 

N c : H s+1 (T c )^H s (r c ), (13) 
see e.g. [14, 19, 27], and they satisfy the classical closed-surface Calderon formula 

N c S c = -^ + K c (14) 

in H S (T C ), where K c is a compact operator. 

While (12) and (13) do not apply to open surfaces, the solutions fi and v of the open-surface 
integral equations (8) and (9) enjoy significant regularity properties. In particular the recent 
results [21] establish that if the arc T and the right-hand-side functions / and g in (8) and (9) are 
infinitely differentiable, we have 

a 

» = - (15) 

to 
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and 

v = P-u, (16) 

where a and j3 are C°° functions throughout T, and where, letting d = d(r) denote any non-negative 
smooth function defined on T which for r in a neighborhood of each end-point equals the Euclidean 
distance from r to the corresponding end-point, oj is a function defined on V such that oj/y/d is 
C°° up to the endpoints. The singular behavior in these solutions is thus fully characterized by the 
factors d 1 / 2 and d~ l l 2 in equations (15) and (16), respectively. 



2.2 Generalized Calderon Formula 

In view of equations (15) and (16), for any non- vanishing function oj(r) > such that ojj\fd is C c 
up to the endpoints of T, we define the weighted operators 



S w [a] = S 



a 

LOJ 



(17) 



and 

N w [/3] = N[/3-u,], (18) 

and we consider the weighted versions 

S w [a] = / (19) 

and 

H U \P] = 9 (20) 

of the integral equations (8) and (9); clearly, in view of the discussion of the previous section, for 
smooth T and smooth right-hand-sides / and g, the solutions a and f3 of (19) and (20) are smooth 
up to the endpoints of T. 

Without loss of generality we use a smooth parametrization r(t) = (x(t),y(t)) of T defined in 
the interval [—1, 1], for which r(i) = is never zero. For definiteness and simplicity, throughout 

the rest of the paper we select oj, as we may, in such a way that 

w(r(t)) = Vl-i 2 - (21) 
The operators S w and N w thus induce the parameter-space operators 

S u [p]{t) = fG k (r(t),r(f)) -?B= r(t')dt>, (22) 

and 

i\UV#) = lim — / G k (r(t) + zn r(t) ,r(t')) ^(t')r(t'Wl - t' 2 dt'; (23) 

2-s>o+ oz j_i on r ((/) 

defined on functions <p and tp of the variable t, — 1 < t < 1; clearly, for (^(t) = a(r(t)) and 
ip(t) = /3(r(t)) we have 

S u [a](r(t)) = S u [<p](t) (24) 

and 

N w t8](r(*)) = JV w [^](t). (25) 
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and 



In order to proceed we further transform our integral operators: using the changes of variables 
t = cos 9 and t' = cos 9' and, defining n# = n r ( cos g) and using (24) and (25), we re-express 
equations (19) and (20) in the forms 

S[0\ = f (26) 

and 

N$] = g, (27) 

where S and N denote the operators 

S[i\{9)= f G k (r{cos9),r{cos9'))-f(9')T{cos9')d9' (28) 
Jo 

N[j}(9)= lim ^- [ ■ 7 ^-G k (r(cos9) + zn e ,r(cos9'))~ / (9 > )T(cos9')sm 2 9 > d9', (29) 

2^0+ OZ Jq OYIqi 

and where 

j{9) = f{v{cos9)) , 0(0)=0(r(cos0)); (30) 
clearly, the solutions of equations (22)-(27) are related by 

<p{0) = f{cos 9) , ${6) = ^(cosfl). (31) 

In view of the symmetries induced by the cos 9 dependence in equations (28) through (30), it 
is natural to study the properties of these operators and equations in appropriate Sobolev spaces 
H^(2tt) of 2-7T periodic and even functions defined below; cf. [8,36]. 

Definition 2 Let s G R. The Sobolev space H^(2tt) is defined as the completion of the space of 
infinitely differentiable 2ir-periodic and even functions with respect to the norm 

oo 

\\vf s = |a | 2 + 2^m 2 > m | 2 , (32) 

m=l 

where a m denotes the m-th cosine coefficient of v: 



1 

v(9) = -ao + ^2 a m cos(m9). (33) 

m=l 

Clearly the set {cos(ra^) : n € N} is a basis of the Hilbert space H^(2tt) for all s. 

The main purpose of this paper is to establish the following theorem. Here and throughout 
this text an operator L between two Hilbert spaces is said to be bicontinuous if and only if L is 
invertible and both L and L _1 are continuous operators. 

Theorem 1 The composition NS defines a bicontinuous operator from H^(2tt) to H^(2n) for all 
s > 0. Further, this operator satisfies a generalized Calderon formula 

NS = J r + K, (34) 

where K : H^(2tt) — > H^(2it) is a compact operator, and where Jq : H^(2tt) — > H^(2n) is a 
bicontinuous operator, independent of k, with point spectrum equal to the union of the discrete set 
Aqo = {Aq = — ^p, X n = — ^ : n > 0} and a certain open set set A s which is bounded away 
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from zero and infinity. The sets A s are nested, they form a decreasing sequence, and they satisfy 
P| S>0 A S = { — j}, where A s denotes the closure of A s . In addition, the operators 

S : H°(2it) -»• F| +1 (2vr) and (35) 

TV: ^ +1 (2vr) H' e (2n) (36) 

are bicontinuous. 

We thus see that, through introduction of the weight uj and use of spaces of even and 2tt periodic 
functions, a picture emerges for the open-surface case that resembles closely the one found for 
closed-surface configurations: the generalized Calderon relation (34) is analogous to the Calderon 
formula (14), and mapping properties in terms of the complete range of Sobolev spaces are recovered 
for S and TV, in a close analogy to the framework embodied by equations (12) and (13). 

In the remainder of this paper we present a proof of Theorem 1. This proof is based on a number 
of elements, the first one of which, presented in Section 3, concerns the operator Jq in (34) — which 
corresponds, in fact, to the zero-frequency /straight- arc version of Theorem 1. 

3 Straight arc at zero frequency: operators Jq and Jq 
3.1 Boundedness of J and link with the Cesaro operator 

In the case in which T is the straight-arc [—1, 1] and k = 0, S reduces to Symm's operator [8, 36] 

S o [0\(9) = -— [ ln|cos#- cos 6'\<p(6)de, (37) 
27r Jo 

whose diagonal property 

f n = Q 

5 [e n ] = X n e n , X n = < J_ (38) 
I 2n ' n — 1 

in the cosine basis {e„ : n > 0} of H^(2tt) (e n (9) = cosnO) is well documented [22]. 

Remark 1 Equation (38) shows that So is a bounded and continuously invertible operator from 
H^(2tt) to H^ +1 (2ir) for all real values of s. 

The corresponding zero-frequency straight-arc version of the operator N, in turn, is given by 

N $\(6) = ^ ]im^2 J \n\{cos 6 - cos 9' f + z 2 \^{e') sin 2 e'dO'. (39) 

Following [13, 19, 26], we express A^o in the form 

A> = D S f (40) 

where 
and 

To[m) = ^m0)sm9); (42) 
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for the sake of completeness, a derivation of the general curved-arc arbitrary-frequency version of 
this relation is provided in Appendix A. 

Note that, in contrast with the closed-arc case [19, p. 117], the expressions (40) through (42) 
contain the vanishing factor sin 9 and the singular factor 1/ sin 9; in particular, it is not immediately 
clear that the operator N maps H^ +1 (2tt) into H^(2tt). To establish this fact, and, further, to 
initiate our general study of the composite operator NS, we start by establishing the following 
result for the straight-arc zero-frequency version 

J = N S (43) 

of the operator NS. 

Lemma 1 The composition Jq = NqSo defines a bounded operator from H^(2tt) into H^(2tt) for 
all s > 0. 

Proof. We first evaluate the action of Jo on the basis {e n : n > 0}. The case n = is straightfor- 
ward: in view of (38) and (40) we have 

MeoW) = (44) 
For n > 0, in turn, expanding (42) we obtain 

To [ e n] (0) = cos 9 cos n9 — n sin n9 sin 9 
cos(n + 1)6* + cos(n - 1)6* cos(n + 1)9 - cos(n - 1)6* ( 45 ) 
= 2 + n 2 

which, for n > 2, in view of (38) yields, upon application of So, 

x ~ r cos(n + l)# cos(n — 1)9 { cos(n + 1)9 cos(n — 1)9\ . . . . 

= W + W + " I W " w ) ■ ( "- 2) (46) 



and, in view of (40)-(41), 



N o [e n }(0) = -cos 9-— - -cos n9, n > 2, (47) 
Zi sin u z 



which can easily be verified to also hold for n = 1. Using (38) and (44) we thus obtain 

Jo[e n ](0) = \ ' fl siring _ cosn6> _ ^ n ( 48 ) 

I cofa(7 4nsin(9 4 ' 

and, therefore, 

J 0lm —M- C -^C lm) + iZ^H jf* W)d9 , (49) 

where the operator C is defined by 

C[e n ](9) = { 1^ ^^ = J (50) 

and where the last term in equation (49) is obtained by collecting the zero-th order terms, and 
explicitly expressing the zero-th order coefficient of (p as an integral. 
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It is not difficult to verify the integral expression 



cmo) = 



fl(7T -0) 

7rsin0 



i r e i r 

a / i f{ u ) du a / <p(v)di 

V JO 7T - V Jg 



(51) 



for C which, in turn, can be expressed in terms of the Cesaro operator 



C[f](x) 



1 f x 

- / f(u)du 
x Jo 



f{xu)du. 



(52) 



As is known [7], C is a bounded operator from L 2 [0, b] into L 2 [0, b] (the space of square-integrable 
functions over [0, b]) for all b > 0. Further, for all / G Cq°[0, b], all m G N, m > 1 and all x G [0, 6] 
we have 



9 m C[/](x) 



dx r 



< 



u 



l f {m) {xu) 



du 



< 



f^\xu) du) =(C[g](x)Y 



(53) 



where, letting g = \f^\, the first inequality follows from equation (52). Integrating this inequality 
with respect to x and taking into account the boundedness of the operator C in I? we obtain in 
particular 

r2n \d m C[f]{x) " 



/ 

Jo 



Q x m 



dx<M\\f^\\l 2[0 ^ 



(54) 



for some constant M. It follows easily from this inequality that C is a continuous operator from 
H™{2it) into H™{2it) for all non-negative integers m. Letting H m (2n) be the space of 2tt periodic 
functions whose derivatives of order k are square integrable in any bounded set of the line for 
all integers k < m (c.f. [19]) we see that C equals the restriction of some continuous operator 
P : H m (2ir) — > H m {2ir): simply define P to equal C on the subspace of even functions, and to 
equal on the space of odd functions. In view of the Sobolev interpolation result (see e.g. [19, 
Theorem 8.13]), P defines a continuous operator from H s (2n) to H s (2ir) for all s > 0, and thus, 
by restriction, 

C maps continuously H^(2tt) ->■ H*(2it) for all s > 0. (55) 

The proof is completed by noting that the last term in the decomposition (49) defines a continuous 
and compact operator from H^(2tt) into itself. □ 



Remark 2 The decomposition (49) resembles superficially the classical closed-surface Calderon 
formula (14), as it expresses the composition Jo = NqSo as the sum of —I /4 and an additional 
operator. As shown in Section 3(3.3), however, the operator C : H^(2tt) — > H^(2ir) which appears 
in (49) is not compact — and thus, formula (49) fails to enable treatment of the operator Jq by means 
of the Fredholm theory. In Sections 3(3.2) through 3(3.4), however, we do show that the operator 
Jo = NoSo is bicontinuous, and, further, we provide a description of its spectrum: the eigenvalues 
of Jq are tightly clustered around — \, in close analogy with the clustering implied by the identity 
term — | in the Calderon formula (14). 

This observation is not limited to the flat-arc zero-frequency case: in Section 5 we further show 
that a certain operator Jq which is closely related to Jq (and which shares its point spectrum) can 
be used as a substitute of the identity operator I in a generalization, valid for open surfaces, of the 
classical Calderon formula. 
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3.2 Invertibility of J 

We proceed to show that the continuous operator Jo admits an inverse which is also bounded from 
H^(2tt) into Hg(2ir). Since the decomposition (40) is not directly invertible on a term by term 
basis (To and Dq are not invertible), we first give an explicit form of the inverse Jq 1 . 

Lemma 2 The operator J : H^(2ir) — > H^(2n) which, according to equations (40) and (43) is 
given by 

J = D S f S , (56) 

is bijective, with inverse 

J^ 1 = -4S^CS f . (57) 

Proof. 

Since the rightmost factor Sq in equation (56) is a diagonal operator, we consider the next operator 
from the right in this product, namely, To. In view of (42) and (50) we clearly have the useful 
identity 

^»]={°;, III. < 58 > 

A corresponding identity for the composition CTq follows from equation (51); in this case we obtain 

Cf = I. (59) 

Having obtained a partial inverse for the factor To in (56) and since the next factor So from the 
right is, once again, a diagonal operator, we consider now the leftmost factor in equation (56): the 
operator Dq. In view of (41) we have 

Do[e n ](e) = { JLi TTl ( 6 °) 

Clearly, this operator bears a simple relationship with the operator C that gives the partial inverse 
of To: in view of (38) and (50) we have 

A> = -ic(V) 2 . (61) 

In sum, to obtain the inverse of Jo we proceed as follows: multiplying Jo on the right by S^C 
we obtain an operator that maps eo to and e n to DoSo[e n ]. Thus, considering (59) and (61), we 
further multiply on the right by — 45oTo and we obtain the operator 

-AJ Q S^CS Q f Q (62) 

which, in view of the fact that the image of ScTo is orthogonal to eo (as it follows easily from 

equations (38) and (45)) maps e n to — 47Jo (<§o) To[e ra ] for all n > 0. But, in view of (61), this 

quantity equals CTo[e n ] which, according to (59), equals e n . In other words, the operator (62) 
maps e n to e n for n = 0, 1, 2 . . . — and, thus, 

Jo = -AS^CSofo (63) 

is a right inverse of Jq. 
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Conversely, since in view of equations (56) and (61) Jq can be expressed in the form 



we have 



Jo — —^CS 1 f So, 



IqJq — S 1 CSqTqCSq 1 TqSq. 



(64) 



(65) 



As noted above, the image of To is orthogonal to en, and thus, since So is a diagonal operator, 
the same is true of the operator S q ~ 1 TqSq. Equation (58) can therefore be used directly to obtain 



f CS 1 f S [e n ] = S 1 f S [e n }, for all n > 0. 
Clearly then, equation (65) can be reduced to 

InJo = S 1 CTqSq, 
and making use of (59), we finally obtain 

-Wo = I, 

as desired. The proof is now complete. □ 

Lemma 3 For all s > 0, the inverse Jq 1 is a bounded mapping from H^(2tt) into H^(2tt). 
Proof. In view of equations (37), (45), (46) and (50) we obtain 



CS f [e n ](9) = { 



l 

2> 
cos 9 
4 ' 

sin(rt+l)0 sin(n— 1)9 f sin(n+l)0 sin(n— 1)9 



4(n+l) 2 sin 8 4(n-l) 2 sin t 



4(n+l) 2 sin 9 4(n-l) 2 sin t 



n = 
n = 1 

n > 2, 



which we re-express as 
where 



CS T [e n ]=F[e n ] + G[e n ], 



FleJ = I 



i 

2' 



n = 
n = 1 

{ cs Q f n+1 V n ' 1 ^ n ^ 2 > 



and where G = CSqTq — F, or, expanding the sine terms in equation (69), 



G[e n ](6) 



1 



0, 



M » COS ( t^jt + t^jjj) +n4^C0M 



(n+1) 2 (n-1) 2 y " sine wo " \Jn+l) 2 (n 

Using (50), the linear operator G can be re-expressed in the form 

G[e n }(9) = R [e n }(9)+ cos 9C [Ri[e„]l (9), n>2 



n = 0,1 
n > 2. 



where the diagonal operators i?o and R\ are defined by 

1 I n n 



Role 



+ 



4 V(n + 1) 2 (n-1 



e„, n > 2 



(66) 
(67) 
(68) 



(69) 

(70) 
(71) 



(72) 



(73) 



(74) 
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and 

Ri[e n ] = 7 — j T\2? T\2 e «' n - 2 ( 75 ) 

(n + l) z {n — \) A 

respectively. Clearly both Rq and R\ are bounded operators from H^(2tt) to H^ +1 (2ir), and, thus, 
in view of (55) and (73), so is G. Invoking once again equation (55) we see from (71) that F is a 
bounded operator from H^(2tt) into -ff| +1 (27r), and in view of (70), the same holds for the operator 
CSqTq. It follows from (57) and (38) that Jq 1 is a bounded operator from H^(2tt) into H^(2tt), as 
claimed. □ 

Corollary 1 For all s > 0, the operator Nq defines a bicontinuous mapping from H* +1 (2t[) to 

Proof. This follows directly from Remark 1, equation (43), and Lemmas 1, 2 and 3. □ 
3.3 Point Spectrum of J 

Having established boundedness and invertibility, we conclude our study of the operator Jq by 
computing its eigenvalues. 

Lemma 4 For any s > 0, the point spectrum a s of Jq : H^(2w) — > H^(2n) can be expressed as the 
union 

a s = A s uA 00 , (76) 

where A^ is the discrete set 

( -111! n = Q 

Aoc = {A n :n = 0,l,...,oo}, A n = _i 4 _' x " (77) 

I 4 4n' n>U ' 

and where A s is the open bounded set 

A s = L = (\ x + i \ y )eC : 4s + 2< { ~fl^l\ 2 \. (78) 



Proof. We start by re-expressing equation (48) as 

f In 2 

M e n}(0) = < si 4 n(n+l)e , cosw g cosn 8 ^ ^ n (79) 
^ An si: 

Then, making use the well-known expansion 

sin(n + 1)0 



[u ' 2 n 

j cosnfl cos nV ~ n 

4risin0 4n 4 ' '* ^ U " 



sm( 



( P 

E (2 - 5o k ) cos 2k0, n = 2p 

k=0 

2 Yl cos(2fc + 1)6, n = 2p+l, 

k=0 



(80) 



(which expressed in terms of Chebyshev polynomials of the first and second kind is given e.g. in 
equation (40) [15, p. 187] and problem 3 in [22, p. 36]) we obtain 



i P_1 
A n e n — ^ 

J [e. 



Ane„ - ± £ (1 - 5 -f)e 2k , n = 2p,p>0 



2n 2 

k=0 ( gl ) 



1 P_1 

A„e n - ^ E e 2fc+i> n = 2p + l, p>0 

fe=0 
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where the diagonal elements A n are defined in equation (77). Clearly, Jo takes the form of an upper- 
triangular (infinite) matrix whose diagonal terms \ n define eigenvalues associated with eigenvectors 
v n , each one of which can be expressed in terms of a finite linear combination of the first n basis 
functions: v n = J2k=o c k e k- ^ n particular, for all n € N, v n 6 i?|[0, 2ir] for all s > 0. This shows 
that the set of diagonal elements defined in equation (77) is indeed contained in a s for all s > 0. 

As is well known, an upper triangular operator in an infinite-dimensional space can have eigen- 
values beyond those represented by diagonal elements. As shown in [7, Th. 2], for instance, the 
point spectrum of the upper-triangular bounded operator 

oo 

c*M(") = EjfeTT (82) 

k=n 

(the adjoint of the discrete Cesaro operator C) is the open disc |A — 1| < 1. As shown in what 
follows, a similar situation arises for our operator Jq. 

To obtain the full point spectrum of the operator Jo let A G C and / = J2h=o fn e n be such that 
<^o[/] = A/. It follows from (81) that the coefficients f n satisfy the relation 



(4-i)/ n -lf ^ = A/ n , n>l, (83) 
v 4 An JJ 2^n + 2k J ' _ ' 
k=i 

along with 



k=i 

Equation (83) is equivalent to 



ly jn±?L = f n{ _ J_ _ 1 A) n>1 (85) 

2 ^ n + 2k J y An 4 ; ' ~ ' v 1 

k=i 

which, by subtraction, gives 

Therefore, the coefficients of / must satisfy 



< 



fn+2= fn 7^T- , n>l 
\ 2 (n+2) / 

OO 

|E#= /o(-^-A), n = 0. 
fc=i 



(87) 



where, in order to simplify the notations, we write 

z = 8A + 2. 

It is clear from equation (87) that the zero-th coefficient is determined by the coefficients of even 
positive orders, and that the sequence f n for n > 1 is entirely determined by f\ and f2- 

Clearly, there are no elements of the point spectrum for which Re(z) > 0, since for such values 
of z the resulting sequence f n is not square summable (that is, 2^ |/n| 2 = oo)- Note that the set 
of vectors {v n } associated with the discrete eigenvalues A n = — \ — in turn, are recovered by 
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setting z 



\. To determine all of the elements of the point spectrum with Re(z) < we study- 



separately the odd and even terms in the sequence (87). We start with the sequence q n = _/2n> 
which satisfies the recurrence relationship 



Qn+i — Qn 



z + 



n>\. 



z — 



(89) 



71+1 



Let z = —x + iy with x > 0, and assume without loss of generality, that q\ = 1. Then 



rt— 1 



n/ k=l 



and it follows that 



ln\q n \ = In 


z - 1 


1 ™ _1 / 

+ E>» 

fc=l V 


= In 


z-l 

n 


1 71—1 

+ EM 

fc=i 



1 ' 



[X 



n > 1, 



i) 2 + 2/ 2 



(90) 



(x + |) 2 + y 2 

1 - r(x,y,k)) 
l + r(x,y,k) 



(91) 



where 



For large k, we have 



and thus 



r(x,y,k) 



2x 



k(x 2 + y* + ±y 



In 



1 - r(x,y,k) 



l + r(x,y,k)J k(x 2 + y 2 
2x 



1. 



ln|g„| = 5 ^lnn + M + 0(-), 

x z + y z n 



where M is a constant. The absolute value of q n is thus asymptotically given by 

1 \ 



O 



nx 2 +y 2 



(92) 

(93) 
(94) 

(95) 



as n — > oo. It follows that, for any s > 0, the set of points (x,y) in the half plane such that the 
sequence ^n 2s \q n \ 2 < oo is exactly defined by the equation 



2s 



4x 



< -1. 



x 2 + y 2 

The analysis for the odd-term sequence p n = fm+\ can be carried out similarly, since 

l 

/ -4 -r - 

/ 7i 

P n+ 1 = P„ 



(96) 



i_ 

n+l+ 



(97) 



which essentially amounts to replacing k by k + \ in equations (91) and (92). The convergence 
condition (96) thus applies to p n as well, and it follows, in view of equation (88), that the set A s 
defined by (78) contains all the eigenvalues of Jq not contained in A M . □ 
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Corollary 2 The operator C : H^(2tt) — > H^(2n) is not compact. 

Proof. This follows from the decomposition (49) of Jo and the fact that Jo admits a spectrum 
that is not discrete. □ 

Remark 3 Using polar coordinates (r,6) around the point (— 4,0) it is easy to check that 



f/> 1 ^, . „ „ cos 6 „ Ttt 37rl ] 

A s = j (Ax + i\ y ) € C : A x + - = rcosfl, A a = rsm0, < r < -— — , 6> € -,y J 



(98) 



Clearly then, for s > s' , A s £ A s / ; and we aave P| S>Q A S = 0, u>/iiZe i/ie intersection of the 
closures is given by f] s>0 A s = {— |}. ^4/so ; /or aZZ s > 0, dist(a s ,0) = — |, and max^ g(Ts |A| < |. 

therefore follows that a s is bounded away from the zero and infinity. In view of Theorem 1 
and Section 3.4, this is a fact of great significance in connection with the numerical solution of 
equations (26) and (27) by means of Krylov- sub space iterative linear- algebra techniques; see [11] 
for details. 

3.4 The operator Jq 

In our proof of Theorem 1 we need to consider not Jo but a closely related operator, namely 

31 = NZ§Z (99) 

where defining (in a manner consistent with equation (106) below) Zq[j](9) = 7(#)r(cos 6), we have 
set 

S^] = S Z [ 1 ], (100) 

and 

N t [ 7 ] = Zq 1 N [ 1 }. (101) 

It is easy to generalize Remark 1, Corollary 1 and Lemmas 1 through 4 to needed corresponding 
results for Sq, Nq and Jq] these are given in the following Theorem. 

Theorem 2 Let s > 0. Then, 

(i) The operator Sq : 7J|(27r) — > iJ| +1 (27r) is bicontinuous. 

(ii) The operator Nq : 2J| +1 (27r) — >■ H^(2ir) is bicontinuous. 
(Hi) The operator Jq : Uf (27r) — > H^(2tt) is bicontinuous. 

(iv) The point spectrum of Jq : H^(2tt) — > H*(2ir) is equal to the point spectrum a s of Jq. 

Proof. In view of (100), (101), the ensuing relation 

J t q = Zq\JqZq, (102) 

and the fact that r is smooth and non-vanishing, the proof of points (i), (ii) and (Hi) is immediate. 
Equation (102) also shows that (A, v) is an eigenvalue-eigenvector pair for Jo if and only if (A, Zq 1 [v]) 
is an eigenvalue-eigenvector pair for Jq, and point (iv) follows as well. □ 
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4 General Properties of the Operators S and N 

The proof of Theorem 1 results from a perturbation argument involving Theorem 2 and the results 
established in this section on the regularity and invertibility of the operators S and TV defined by 
equations (28) and (29). 

4.1 Bicontinuity of S 

We seek to show that for all s > the operator S defined in equation (28) is a bicontinuous mapping 
between H*(2ir) into -ff| +1 (27r). This is done in Lemmas 5 and 7 below. 

Lemma 5 Let s > 0. Then S defines a bounded mapping from H^(2tt) into H^ +1 (2n). Further, the 
difference S — Sq (see equation (100) ) defines a continuous mapping from H*(2tt) into iJ| +3 (27r). 

Proof. In view of equation (4) and the expression 

Hl{z) = -J (z) \tv(z) + R{z) 

7T 

for the Hankel function in terms of the Bessel function Jo(z), the logarithmic function and a certain 
entire function R, the kernel of the operator 5 W (equation (22)) can be cast in the form 

G k (r(t), r(0) = A 1 {t, t') In \t - t'\ + A 2 (t, t>), (103) 

where A\(t, if) and A 2 (t, if) are smooth functions. Further, since Jq(z) is given by a series in powers 
of z 2 , it follows that for all m € N, the function A\ can be expressed in the form 

^ TO+3 

Ai(t,0 = + ^(W - tT + (* - tT + ^m+3(t,t'), 

n=2 

where A m+ s(t,t') is a smooth function of t and t'. The operator S in equation (28) can thus be 
expressed in the form 

S[0\{6) = f%[0\(9) + ^a„(cosfl) / (cosfl' - cos#) n ln| cos8 - cos 8' \p (8' )t (cos 8' )d8' 

n=2 J ° (104) 



+ 



/•7T 

/ A 3 (cos 8, cos 8')<p{8')t(cos 8')d8', 
Jo 



where A3 (cos 8, cos #'), which contains a logarithmic factor, belongs to C m+3 ([0, 2-zr] x [0, 27r]). 

Clearly, for n > 2, the second derivative d 2 /d8 2 of the product (cos 8' — cos 8) n In | cos 8 — cos 8'\ 
can be expressed as a product Pi (cos 8, cos #') In | cos 8 — cos 8' \ + i-2(cos 0, cos 8') where P\(t, t') and 
P2(t,t') are polynomials. Collecting terms with the common factor cos 1 8' we then obtain 

,2 m+1 -Tj- 

— 2 (-S-Sj) = be(.cos8)S Z e [ip](8)+ / ^(cos^cosfl'^'Mcosfl')^', (105) 

e=o 

where be(cos8) is an even smooth function, where the operator Zg : H^(2tt) — > H^(2tt) (s G M) is 
given by 

Z*[7](0') = cosV r(cos0') 7(0'), (106) 
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and where A 4 (cos 9, cos 9') £ C m+1 ([0,27r] x [0,2vr]). Now, in view of Remark 1, the first term on 
the right-hand-side of equation (105) defines a bounded operator from H^(2n) into H^ +1 (2tt). On 
the other hand, the derivatives of orders k < (m + 1) of the second term on the right-hand-side 
of (105), all reduce to integral operators with bounded kernels, and thus map L 2 [0, 2ir] continuously 
into L 2 [0, 2ir]. It follows that the second term itself maps continuously H®(2it) (and hence H™(2ir)) 
into H™ +1 (2ir), and the lemma follows for integer values s = m. The extension for real values s > 
follows directly by interpolation [19, Theorem 8.13]. □ 

The following lemma and its corollary provide a direct link between the spaces H*(2ir) under 
consideration here, and the original space H~z(T) appearing in equations (10). 

Lemma 6 Let s > 0, and assume (p 6 H^(2tt). Then the function 

«,(£) = - {p{e)e-* cose de. (107) 

satisfies 

"" K)|2 T dt < oo. (108) 



/ 



(i + Kl 2 )* 

Proof. Using the L 2 [0, 7r] -convergent cosine expansion 

oo 

<f(9) = ^2a n cos9 (109) 

n=0 

we obtain 

W (£\ = Y— cosn9e-^ cos8 d9. (110) 

Since 



Jo 



cos n9e-'^ cos 6 'al9 = \ T e 1 ™ 6 e'*™ 56 dO = \e^ f* e ~ ine e'^ 9 d9 = m n J n (-0, (111) 



(where, denoting by J n {€) the Bessel function of order n, the last identity follows from [16, 8.411 
p. 902]), we see that equation (110) can be re-expressed in the form 

w(0 = £ i n a n J n {-i) = f; (Vl + n^ i n a n ) ( %t^ \ . (112) 

In view of the Cauchy-Schwartz inequality we thus obtain 

koi 2 < (Ea+^Ki 2 ) (E t0) * (e^ + woi 2 ) \ml (H3) 

Since < |f |/(1 + l^ 2 ) 1 / 2 < 1, it follows that 



:i + iei 2 )^ V&V n2s ^(i 



^(Ol 2 M \ , /• |J (OI 



rde + ' ' i <£ \\<p 



~ II 2 



en j--^ v^tv n ^(i + iei 2 )^ / ^(i + ici 5 

Jffl /" kn(Q| 2 ,A . /• |Jo(Q| 2 ,A ,,,-,,2 



(114) 
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Further, in view of [16, 6.574, eq 2.], the integral involving J n can be computed exactly for n > 1: 



/ 

JR 



d£ = -. (115) 



lei s n 



It thus follows that 



where 



/ 

JR 



\w(0\ 



2 



(i + iei 2 ) 5 



-d£ < C s \\(p\\i < oo (116) 



n=i n •Mi + iei 2 ) 2 



□ 



Corollary 3 Let s > 0, <p e H*(2ir), ip(t) = <p (arccos(i)) , <^ : [-1, 1] ->■ C, a(p) = ^(r _1 (p)) and 
W(p) = w(r _1 (p)). Then, the function F = ^ is an element of H~z(T). 

Proof. It suffices to take show that / = ip/u G H~? [— l i 1] for the case T = [—1, 1]. Extending / 
by outside the interval [—1, 1], the Fourier transform of / is given by 



7(0 = f°° f(t)e-^dt= [ y?(t)6 " 1 dt= r \p(9)e-^ cose d9, (118) 

J-OD J-l UJ { t ) JO 

since oj(t) = y/l — t 2 in the present case. The Corollary now follows from Lemma 6. □ 

Lemma 7 For all s > the operator S : H^(2ir) — > H^ +1 (2tt) is invertible, and the inverse 
S^ 1 : Hl +1 (2n) — > H^(2tt) is a bounded operator. 

Proof. Let s > be given. From equation (38) it follows that S : H*(2ir) ->■ H* +1 (2ir) is a 
continuously invertible operator. The same clearly holds for Sq as well, and we may write 

S = S^I+^y\s-S^y (119) 

It follows from Lemma 5 that the operator (5J)~ 1 (,§ — Sq ) is bounded from H^(2tt) into H^ +1 (2tt), 
and therefore, in view of the Sobolev embedding theorem it defines a compact mapping from H^(2tt) 
into itself. Further, in view of Corollary 3 and the injectivity of the mapping (10) it follows that 
the operator S : H*(2ir) — > H* +1 (2ir) is injective, and therefore, so is 



e 

(S^ 1 S = I+ (S^y 1 (S - S5) : ff|(27r) ->• H' e (2n). (120) 

A direct application of the Fredholm theory thus shows that the operator (120) is continuously 
invertible, and the lemma follows. □ 
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4.2 Bicontinuity of iV 

To study the mapping properties of the operator N we rely on Lemma 8 below where, as in [26], 
the operator N is re-cast in terms of an expression which involves tangential differential operators 
(cf. also [13, Th. 2.23] for the corresponding result for closed surfaces). The needed relationships 
between normal vectors, tangent vectors and parametrizations used are laid down in the following 
definition. 

Definition 3 For a given (continuous) selection of the normal vector n = n(r) on T, the tangent 
vector t(r) is the unit vector that results from a 90° clockwise rotation of n(r). Throughout this 
paper it is further assumed that the parametrization r = r(t) of the curve T has been selected in 
such a way that 

t(r(t)). (121) 

Lemma 8 For tp £ C°°(F), and fort £ (—1, 1), the quantity iV w [ip](i) defined by equation (23) can 
be expressed in the form 

N u [<p]{t) = N*[<p]{t) + N?[<p](t) (122) 

where 

= k 2 J G k {v{t),v{t')) tptf) r(t') Vl-t' 2 n t • u t , dt', (123) 

and where 

K v m) = ^| (/' G fc (r(t),r(f )) (W) VT^) dt'^j . (124) 

Proof. See Appendix A. □ 

In order to continue with our treatment of the operator N we note that, using the changes of 
variables t = cos 9 and t' = cos 6' in equations (123) and (124) together with the notation (31), for 
Lp G C°°(r) and for 9 G (0,vr) we obtain 

N[0\ = N 9 &] + N pv &], (125) 

where 

N 9 [0\{9) = k 2 [ G k (r(cos9),r(cos9')) (p[&) r{cos9') sin 2 9' n e ■ n e , d9', (126) 
Jo 

and where, taking into account equations (41) and (42), 

^ PV ^ = 7(^) {°oSn) [0\{O), (127) 

with 

* 9) = « W1(9) - (128) 

Lemma 9 Let s > 0. The operator N pv defines a bounded mapping from H^ +1 (2tt) to H^(2tt). 
Further, the difference (N pv — Nq) (see equation (101) ) defines a bounded mapping from H* +1 (2ir) 
into H s + 1 (2t:). 



dr. . 



dr 

dt 
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Proof. Using (40) (101) and (127) we obtain 

N*°m = Nsm + ^^Dois - smm. (i 29 ) 

As shown in Theorem 2 the operator Nq : Hl +l (2ir) — > H^(2tt) on the right-hand side of this 
equation is bounded. To establish the continuity of the second term on the right-hand side of 
equation (129) we first note that, in view of equation (42), the operator To : H^ +1 (2n) — > H^(2n) 
is bounded, and therefore, so is Tq~. Further, as shown in Lemma 5, the operator (S — Sq) maps 
continuously H^(2n) into -fT| +3 (27r) so that, to complete the proof, it suffices to show that the 
operator D maps continuously i7| +3 (27r) into H^ +1 (2ti). But, for tp <E H* +3 (27r) (s > 0) we can 
write 

and since the zero-th order term in the cosine expansion of Jp-^ vanishes, in view of (50) we have 

d 2 j>~ 

It therefore follows from (55) that the second term in (129) is a continuous map from H^ +1 (2ir) 
into H e s (27r), that is, {N pv - A> T ), as claimed. □ 

Corollary 4 For all s > the operator N can be extended as a continuous linear map from 
H^ +1 (2tt) to H^(2tt). Further, the difference N — Nq defines a continuous operator from H^ +1 (2n) 
to H s e +1 {2n). 

Proof. From equation (126) we see that N 9 has the same mapping properties as S (Lemma 5), 
namely 

N 9 : H s e {2ir) -> H s e +1 {2n) is continuous. (130) 
In view of Lemma 9 it therefore follows that the right hand side of equation (125), 

N 9 + N pv : H s e +1 (2ir) -»• # e s (2vr), (131) 

is a bounded operator for all s > 0. Equation (125) was established for functions (p of the form (31) 
with tp e C°°(r). But the set of such functions tp is dense in H^ +1 (2tt) for all s > — as can 
be seen by considering, e.g., that the Chebyshev polynomials span a dense set in H s+1 [— 1, 1]. It 
follows that N can be uniquely extended to a continuous operator from i?| +1 (27r) to H^{2ir), as 
claimed. Finally, N - = N 9 + (N pv - Nq) is continuous from H s e +1 (2Tr) into J ff| +1 (2vr), in view 
of equation (130) and Lemma 9. 
□ 

The following lemma establishes a link between the domain of the unweighted hypersingular 
operator N considered in [32] (equation (11) above) and the corresponding possible domains of the 
weighted operator iV (equation (36)); cf. also Corollary 3 where the corresponding result for the 
domains of the operators S and S is given. 

Lemma 10 Let tj) belong to J ff| +1 (27r) for s > 0, ip(t) = ^(arccost), tp : [-1,1] ->■ C, /3(p) = 
ip (r _1 (p)) ; W(p) = oj (r _1 (p)). Then the function G = W(3 is an element of H^{Y). 



D $] = C 
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/ 

JR 



Proof. It suffices to show that g = uij} G H 1 * [— 1, 1] for the case T = [—1,1]. Extending g by 
outside the interval [—1,1], the Fourier transform of g is given by 

9(0 = J ^{ty-^oo^dt = tf;(cos 9)e~^ cose sin 2 9d9, (132) 
since oj(t) = \/l — t 2 in the present case. Integrating by parts we obtain 

9(0 = -^Mco S e)sm0}e-^ cose de. (133) 

It is easy to check that ^{V>(cos 9) sin#} = -§g{ip(9) sin 9} is an element of H*(2n) and, thus, in 
view of equation (133) together with Lemma 6 we obtain 

\q(£)\ 2 £ 2 

m?jl < oo (134) 

It thus follows that the second term on the right-hand side of the identity 

/ \9(0?(l + \^di= I ^-rdi+ [ (135) 
Jr (1 + £ 2 )5 (1 + £ 2 ) 2 

is finite. The first term is also finite, as can be seen by applying Lemma 6 directly to equation (132). 
The function g thus belongs to Hi[—1, 1], and the proof is complete. □ 

Lemma 11 For all s > the operator N : H^ +1 (2tt) — > H^(2ir) is invertible, and the inverse 
N^ 1 : H*(2tt) -)• H^ +1 (2tt) is a bounded operator. 

Proof. In view of Theorem 2, the operator Nq : H'1 +1 {2tt) — > H^(2tt) is bicontinuous, and we 
may thus write 

TV = N5 (l+ (N - N5)\ . (136) 

Since, by Corollary 4, the difference N — Nq defines a bounded mapping from H^ +1 (2tt) into 

H* +1 (2tt), it follows that the operator (iV T ) 1 (iV-iVg") is bounded from # e s+1 (2vr) into # e s + 2 (2vr) 

and, in view of the Sobolev embedding theorem, it is also compact from H^ +1 (2ir) into H^ +l (2ir). 
The Fredholm theory can thus be applied to the operator 



1 + 



(N5)\N-N5). (137) 



This operator is also injective, in view of Lemma 10 and the bicontinuity of the map N in equa- 
tion (11), and it is therefore invertible. The Lemma then follows from the bicontinuity of the 
operator of Nq. □ 



5 Generalized Calderon Formula: Proof of Theorem 1 

Collecting results presented in previous sections we can now present a proof of Theorem 1. 
Proof. The bicontinuity of the operators S, N and NS follow directly from Lemmas 5, 7, 11 and 
Corollary 4. To establish equation (34), on the other hand, we write 

NS = N^S^ + K = J t + K (138) 
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where as shown in Theorem 2, is bicontinuous, and where 

k = N(S-S^ + (N-N^)S^. (139) 

In view of Lemma 5, Corollary 4 and Theorem 2, the operator K maps H^(2tt) into H^ +1 (2tt) and 
is therefore compact from H^(2w) into H^(2ir). The proof is now complete. □ 

A Proof of Lemma 8 

Proof. Assuming ip G C°°(T), we define the weighted double layer potential by 

T> u [a](r) = ! dG } r,r 'K (r')oj(r')de', r outside V, (140) 
Jr on r < 

which, following the related closed-surface calculation presented in [13, Theorem 2.23], we rewrite 
as 

D w [a](r) = -div r E[a](r) (141) 

where 

E[a](r) = J^G k {v,v')a{v')u{v')rv{v')dl'. (142) 

Since E[a] = E = (.E x , E y ) satisfies AE + /c 2 E = 0, the two dimensional gradient of its divergence 
can be expressed in the form 

/ Q \ 

grad divE = -A; 2 E + —curl E, -—curl E , (143) 
\dy dx J 



where the scalar rotational of a two-dimensional vector field A = (A x ,A y ) is defined by curl A = 
(l^T ~ TJtfO- Since curlr ( n ( r ')G'(r,r / )) equals -t(r') • V r /G(r,r') (see definition 3), we obtain 

curl r E[a](r) = - J* dG ^W a{v(t'))u>{v{t'))dt' . (144) 

Therefore, taking the gradient of (141), letting f(t') = a(r(t')), using (21), integrating (144) by 
parts, and noting that the boundary terms vanish identically (since \/l — t' 2 = for t! = ±1), we 
see that 

gradD w [«](r) = fc 2 ^G(r,r')«(r')n(r / ^(r / K / - , (145) 

A(r) = j^G(r,r(t'))^(v(t')VT^)dt'. (146) 

In view of the continuity of the tangential derivatives of single layer potentials across the integration 
surface (e.g. [13, Theorem 2.17]), in the limit as r — > r(t) 6 V we obtain 

dA(r(t)) dA(r(t))\ _ 1 dA(r(t)) ^ (M?) 



where 



dy ' dx J ' T (f^ dt ' 

and the decomposition (122) results. □ 
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B Asymptotic behavior of NS[1] 



In this section we demonstrate the poor quality of the composition NS of the unweighted hyper- 
singular and single-layer operators by means of an example: we consider the flat arc [—1,1] at zero 
frequency (NS = NoSo). In detail, in Section B.l we show that the image of S is not contained 
in the domain of N (and, thus, the formulation NS cannot be placed in the functional frame- 
work [32,33,35]), and in Section B.2 we study the edge asymptotics of the function NS[1] which 
show, in particular, that the function 1, (which itself lies in H s [— 1, 1] for arbitrarily large values 
of s) is mapped by the operator NS into a function which does not belong to the Sobolev space 
H~z[— 1,1], and, thus, to any space H s [— 1,1] with s > —1/2. 

We thus consider the unweighted single-layer and hypersingular operators which, in the present 
flat-arc, zero- frequency case take particularly simple forms. In view of (6), the parameter- space form 
of the unweighted single-layer operator (which is defined in a manner analogous to that inherent 
in equation (22) and related text) is given by 

i r 1 

S [<p](x) = -— J \n\x-8\ip(8)ds. (148) 

With regards to the parameter-space form Nq of the hypersingular operator (7) we note, with 
reference to that equation, that in the present zero-frequency flat-arc case we have r = (x,0), 
zn r = (0,z) and —d/dn r / = d/dz. Since, additionally, the single layer potential (2) yields a 
solution of the Laplace equation in the variables (x, z), we have 

d 2 r 1 

4irN [ip]{x) = -lim— ^ / <p{s) ln((x - s) 2 + z 2 )ds, (149) 
z^o dx L J_i 

or equivalently, 

AT„M(*) = 1. to i. £ H(x _ s? + z i )is . (150) 

Note that, in view of the classical regularity theory for the Laplace equation, letting z tend to zero 
for —1 < x < 1 in equation of (149) we also obtain, for smooth ip, 

N [<p](x) = -^SoM(x), -Kx<l. (151) 



B.l The operator S 

Integrating (148) by parts we obtain 

-2ttS [1](x) = f d ( S ~ ^ In \s - x\ds = (1 - x) ln(l - x) + (1 + x) ln(l + x) - 2, (152) 
J-i ds 

and therefore 

So[l](x) = ^-(2 - (1 - x)ln(l - x) - (1 + x) ln(l + x)). (153) 

Incidentally, this expression shows that the unweighted single-layer operator does not map C°° 
functions into C°° functions up to the edge; a more general version of this result is given in [34, p. 
182]. 

The following two lemmas provide details on certain mapping properties of the operator So . 
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Lemma 12 The image Sq[1] of the constant function 1 by the operator (148) is an element of 
H*[-l,l]. 



Proof. Let Fi be a closed, smooth curve which includes the segment [—1, 1]. Clearly, the function 

belongs to L 2 (Fi) and therefore to H~2(Ti), so that, according to Definition 1, the constant 1 is 
in the space H~z [— 1, 1]. In view of equation (10), it follows that Sq[1] G H*[—l, 1]. □ 

Lemma 13 The image So\l] of the constant function 1 by the operator (148) is not an element of 

hH-i,i}. 

Proof. In view of (153) and the fact that So[l](x) is an even function of x, integration by parts 
yields 

C e-*%[l](s)dz = ^S [l](l) + -Kt C e-^\n l —^dx. (155) 
Taking into account the identities [16, eq. 4.381, p 577] 

f fthixcos&dx = -\ [si(0 + f] , ^ >Q (156) 

\ fo^ixsin^xdx = -| [C + ln£ - ci(£)] , 

where C is the Euler constant, and where si(£) and ci(^) are the sine and cosine integrals respec- 
tively, (both of which are bounded functions of £ as |£| tends to infinity), it is easily verified that 
the second term in (155) behaves asymptotically as ^jP- as £ tends to infinity. Clearly, the first 
term of (155) decays as 0(|), and therefore 



e~^ x (S [l](x))dx 



-i 



Equation (157) tells us that the function ip : R — > R which equals So[l](x) for x in the interval 
[— 1, 1] and equals zero in the complement of this interval, does not belong to H2 (R), and, thus, 
5 [1] H*[-l, 1], as claimed. □ 

Remark 4 Lemmas 12 and 13 demonstrate that, as pointed out in Section 1, the formulation NS 
of the open-curve boundary-value problems under consideration cannot be placed in the functional 
framework put forth in [32, 33, 35] and embodied by equations (10), (11) and definition 1: the image 
of the operator S is not contained in the domain of definition of the operator N; see equations (10) 
and (11). 



B.2 The combination N S 

While, as pointed out in the previous section, Sq[1] does not belong to the domain of definition of 
Nq (as set up by the formulation (10), (11)), the quantity NqSo[1](x) can be evaluated pointwise 
for I a; I < 1, and it is instructive to study its asymptotics as x — > ±1. 
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Lemma 14 7V So[l] can be expressed in the form 



N So[l](x) = ^ 1 + C(x), (158) 

7T Z (1 — X z ) 

where C G L 2 [-l,l]. 

Proof. In view of (153) we have 

JV 5 [l](x) = -iV [l](x) - ^-iV [ 5 ](x), (159) 

where 

5 (x) = (1 - x ) l n (l - x ) + (1 + x ) ln(l + x). (160) 
For the first term on the right-hand side of this equation we obtain from (151) and (153) 

N [l](x) = - 1 (161) 

7T(1 — J 

To evaluate the second term iVo[g] in equation (159), in turn, we first integrate by parts equa- 
tion (150) and take limit as z — > and thus obtain 



(162) 



No[9]{x) = h{t ln|x - s|ff(s) ^) - Tx L ln|x - s ^ 9{s)ds ) 

= l] VTx ( ln (iTi) ) + ln |x " s| 111 ( Tt£) ds ' 

or 

w °t*> = ^ - /_! (itt) i^*- < 163 > 

Clearly, to complete the proof it suffices to establish that the functions 

C + (x)=p.v. I' Hl ~ s) ds and C-(x)=p.v. C ln(1 + §) ds (164) 
7-1 s-x J-i s-x 

are elements of L 2 [— 1, 1]. 

Let us consider the function £ + for x > first. Re-expressing C + (x) as the sum of the integrals 
over the interval [x — (1 — x),x + (1 — x)] = [2x — 1, 1] (which is symmetric respect to x plus the 
integral over [— 1, 2x — 1] and using a simple change of variables we obtain 



„. , , r~ x ln(l -x-u)- ln(l - x + u) , /■ i+x ln(l - x + u) , , N 

£+ (x )= / — -du + / — -du. 165) 

Jo u J x _ x u 

Letting z = 1 — x and v = |, we see that the first integral in (165) is a constant function of x: 
f z ln(z - u) - \n(z + u) , Z" 1 ln(l - v) - ln(l + u) , 

/ — i '- i Ldu = / — i '- i '-dv = const. (166) 

Jo u Jo v 
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For the second integral in (165), on the other hand, we write 



Mi + -*),», _ ki + A)^ + ln(1 _ x) p * 



/■it: Ml+^) ^ + ln(1 _ x)ln (l±^ 



= Ml±^ + l ^±4dv + ln(l - x) In ( i±£ 



(167) 



1 f in' f ^- ) - In' 2 ) + P= ^±^, + ln(l - x) In (\±± 
2\ \l-xj J J, v(l + v) V ; V 1 -^ 



Since the second term on the last line of equation (167) is bounded for < x < 1, it follows that, 
in this interval, the function C + {x) equals a bounded function plus a sum of logarithmic terms 
and is thus an element of L 2 [0, 1]. Using a similar calculation it is easily shown that C + (x) is 
bounded for —1 < x < 0, and it thus follows that C + G L 2 [—l, 1], as desired. Analogously, we have 
Cr e L 2 [—l, 1], and the lemma follows. □ 

Corollary 5 Let T = [—1,1]. Then NS[1] does not belong to the codomain H^^[— 1,1] of the 
operator N in equation (11). 

Proof. In view of Lemma 14 it suffices to show that the function h{x) = jz^tz does not belong to 
H~^[— 1,1], or, equivalently, that the primitive k(x) = — \ ln^| of h does not belong to H^[—l,l]. 
Clearly, to establish that k g" H^[— 1,1] it suffices to show that the function £(x) = p(x)ln(x) is 
not an element of H2 [0, oof, where p is a smooth auxiliary function defined for x > which equals 
1 in the interval [0, 1] and which vanishes outside the interval [0, 2]. 
To do this we appeal to the criterion [20, p. 54] 

^ e ff2(0,oo) ^ ^GL 2 (0,oo) and/ t~ 2 dt \i(x + t) - £(x)\ 2 dx < 00. 

Jo Jo 

To complete the proof of the lemma it thus suffices to show that the integral 

1=1 t- 2 dt |ln(x + t) -ln(x)\ 2 dx 
Jo Jo 

is infinite. But, using the change of variables u = j we obtain 
and the lemma follows. □ 
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